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ABSTRACT 

Boxy/pcanut bulges are believed to originate from galaetie dises through secular pro- 
cesses. A little explored question is how this evolution would be modified if the initial 
disc was assembled around a preexisting classical bulge. Previously we showed that a 
low- mass initial classical bulge (ICB), as might have been present in Milky Way-like 
galaxies, can spin up significantly by gaining angular momentum from a bar formed 
through disc instability. Here we investigate how the disc instability and the kinemat- 
ics of the final boxy/peanut (BP) bulge depend on the angular momentum of such a 
low- mass ICB. We show that a strong bar forms and transfers angular momentum to 
the ICB in all our models. However, rotation in the ICB limits the emission of the 
bar's angular momentum, which in turn changes the size and growth of the bar, and 
of the BP bulge formed from the disc. 

The final BP bulge in these models is a superposition of the BP bulge formed via 
the buckling instability and the spun- up ICB. We find that the long-term kinematics of 
the composite BP bulges in our simulations is independent of the rotation of the ICB, 
and is always described by cylindrical rotation. However, as a result of the co-evolution 
between bulge and bar, deviations from cylindrical rotation are seen during the early 
phases of secular evolution, and may correspond to similar deviations observed in 
some bulges. We provide a simple criterion to quantify deviations from pure cylindrical 
rotation, apply it to all our model bulges, and also illustrate its use for two galaxies: 
NGC7332 and NGC4570. 

Key words: galaxies: bulges - galaxies: structure - galaxies: kinematics and dynamics 
- galaxies: spiral - galaxies: evolution 



1 INTRODUCTION 

The formation of realistic disc galaxies and their evolution 
has remained a challenging problem to the cold dark matter 
(CDM) paradigm of galaxy formation (|White fc Rccs lOZl). 
Nevertheless, significant development has occurred in this 
context in recent years. Cosmological hydrodynamical sim- 
ulations which include feedback and/or smooth accretion 
of cold gas through cosmic filaments show that an expo- 
nential disc could have assembled either around a classical 
bulge formed through mergers or as a bulgeless disc galaxy, 
and survived through the hierarchical assembly 
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Governato et al.ll2007l : lAgertz et alll201ll : 
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Soon afte r their formation, t he clumpy discs observed at 
high-redshift jCenzel et a il l2006l ) are believed to go through 
rapid dynamical evolution where minor mergers, streams of 
inflowing gas triggering burst of star formation and tur- 



E-mail:saha@mpe. mpg.de 



bulence played dominant roles. As th ese violent processes 
become less imp ortant and frequent (|Hopkins et al.l l2010l : 
iLotz et al]|201ll ). secular processes which operate on longer 
time scales (typically ~ 10s of dynamical time scales) are 
thought to play a signi ficant role in the subsequent evo- 
lution of the galaxies (|Kormendv fc Kennicuttl l2004h . In 
the secular phase, the most efficient way a disc galaxy 
evolves is through forming a bar via disc instability which 
facilitates the redistribution of energy and angular mo- 
mentum between the disc, d ark matt er halo and classi- 
cal bulge jPebattista fc Sellw ood 2000j; lAthanassouiall2003l : 
ISaha et al.ll2012l ). As the bar grows stronger, it goes through 
buckling instability and form boxy /peanut bulges (hereafter, 
BP) as demonstrated in nume r ous N-body simulation stud- 
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m nearly 50% of edge-on disc galaxies (jLiitticke et al.|[2000l ) 
including our Milky Way and the line-of-sight (LOS) stel- 
lar and gas kinematics in many such barred edge-on galax- 
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ies show cylindrical rotati on fKormcndv & IllinKworth"l98^; 
[Bureau fc FreemanI 1 19991 : [Falcon-Ba rroso ct al- .2006a ), as 
has also been noted in many N-body simulations of barred 
galaxies (ICombes et al.]|l990l: ISellwood fc \yilkinsonl Il993l : 



lAthanassoula fc MisiriotisI |2002| : ISaha et alj \2Qli ). Stellar 
kinematics from the recently completed Bulge Radial Veloc- 
ity Assay (BRAVA) surv ey confirms the cylindrical rotation 
in th e Galactic Bulge (jHoward et alJ |2009| : iKunder et alJ 
|2012| ). Recent observational analyses, however, show devia- 
tions from cylindrical rotation in BP bulges, e.g., NGC5746, 
NGC1381 (| Williams et all 120111 ). Deviations from cyhndri- 
cal rotation is also reported in BP bulges of simulated 
galaxies due to variation in the viewing angles of t he bar 
llCombes et al ]|l990l : lAthanassoula fc Misirioti3l2002l ') . How- 
ever, it is not clear whether the projection effect is the cor- 
rect interpretation, or whether a physical mechanism is re- 
quired to explain the deviations from cylindrical rotation 
observed in BP bulges. 

BP bulges are not only formed out of pure axisym- 
metric discs but they can also for m in axisymmetri c discs 
assembled around classical bulges (|Saha et al.ll20l3 ). Such 
classical bulges are presumably have formed as a result 
of vi olent processe s e.g., dissipative collapse (|Eggen et al.l 
Il962f ) or mergers (|Baugh et al.|[l99^ V In observation, one 
does find both barred and unbar red galaxies with classi- 
cal bulges (|Laurikainen et aLll2007l ). Classical bulges are in- 
deed abundant in disc galaxies and they are known to rotate 
l|Kormendv fc Illingworthll 19821: ICappellari et alll2007l '). In a 



previous paper, ISaha et al.1 (|2012l ) showed how an initially 
non-rotating low-mass classical bulge spun-up during the 
secular evolution and discussed the possible implication it 
might have on the final BP bulge. Such a low-mass classi- 
cal bulge [r^ 8% of disc mass) might be present in the Milky 
Way as suggested by the A^-body modelling of BRAVA stellar 
kinematics; models with increased bulge masses are show n 
to produce larger deviation from the data (|Shen et al.ll2O10l ). 

The goal of this paper is to understand the effect of 
angular momentum of a low-mass ICB on the disc instabil- 
ity and kinematics of the final BP bulge. In particular, does 
the cylindrical rotation in BP bulges depend on the angular 
momentum of the ICB and could this possibly explain the 
deviation from cylindrical rotation in observed BP bulges? 
Here, we study the formation and growth of a bar and BP 
bulge in galaxies with a low mass rotating ICB using high 
resolution N-hoAy simulations. We find that the kinemat- 
ics of the final BP bulge varies with time but the long-term 
asymptotic behaviour of the resulting BP bulge does not 
depend on the initial rotation. We show how cylindrical ro- 
tation in BP bulges changes during the secular evolution and 
devise a method to quantify it. 

The paper is organized as follows. Section [2] outlines the 
initial galaxy model and set up for the A'^body simulation. 
Bar growth and the formation of BP bulges are described 
in section |3] Kinematics of boxy bulges are described in sec- 
tion 13] Conclusions are presented in section [S] 



2 GALAXY MODELS WITH INITIALLY 
ROTATING CLASSICAL BULGES 

Equilibrium models of galaxy with I CBs are constructed 
using the self-consistent method of iKuiiken fc Dubinskil 




Figure 1. Radial profile of ToomrG Q for the initial stellar disc. 



0.6 




t (Gyr) 

Figure 2. Top panel showing the time evolution of the bar 
strength in galaxy models with a preexisting spinning classieal 
bulge. Bottom panel shows the net gain of angular momentum 
by the preexisting classical bulge. L^^ (0) denotes the disc specific 
angular momentum at T=0. Black lines in both panels denote 
the model with non-rotating ICB (RCG004-0). Red (RCG004- 
A), green (RCG004-B) and blue {RCG004-C) Unes in both panels 
indicate models with rotating ICBs. 



(|l995h . We construct here three galaxy models with rotat- 
ing I CBs and one with non-rotating ICB RCG004-0 taken 
from lSaha (|2012D for comparison. Each of these initial 

models consists of a live disc, halo and bulge. The initial disc 
has an exponentially declining surface density with a scale- 
length Rd, mass Md and Toomre Q = 1.4 at 2.5Rd (see 
Fig. [1] for the radial variation). The live dark matter halo 
is modelled with a lowered Evans model and the ICB with 
a King model and gen eral details about the model param - 
eters can be found in (|Saha et al.1 l2010l : ISaha et al.ll2012l ) . 
The initial parameters for disc, halo and ICB are identical 
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to model RCG004 of ISaha et al] (|2012h and the circular ve- 
locity curve is the same as presented in Fig. 1 of I S aha et alj 
( 20121). The ICBs are st rongly flattened by the disc potential 
( Barnes fc WhitdllQsA . The ICE in the constructed models 
has a total mass Mb = 0.066Md and the initial ellipticity in 
edge-on projection is et — 0.465. 

The rotating flattened models of ICBs are constructed 
by reversing the velocities of particles with negative angular 
momenta, which rem ains a valid solution of the coUisionless 
Boltzmann equation (|Lvnden-Belllll96^ ) . In this way, three 
models with rotating ICBs which we denote as RCG004- 
A, RCG004-B and RCG004-C are constructed. Their initial 
velocity contours are shown in section [l] Model RCG004- 
A were constructed by randomly changing the corotating 
fraction (i.e., a fraction of bulge orbits corotating with the 
disc) to 0.75 (note, corotating fraction = 0.5 generates a 
non-rotating bulge) and it shows some amount of cylindrical 
rotatio n in the centra l region, although the ICB is axisym- 
metric (|Rowlevl[l98l ). Models RCG004-B and RCG004-C 
were constructed sl i ghtly differently based on the work of 
iDehnen fc Gerhard! (|l993l ). in that we reversed the veloci- 
ties of bulge stars belonging to a certain range of angular 
moment a, especially avoiding the m eridional loop orbits, see 
Fig. 3 of lDehnen fc Gerhardl (| 19931 ). 

The initial {Vm/o')* values for all the model bulges are 
shown in Table [T] where Kn is the maximum velocity and a 
is the average velocity dispersion inside the bulge half mass 
radius. The * sign denotes that values are nor malized to th e 
corresponding oblate isotropic rotator model (|Binnevlll97"i ). 

We scale the models such that Md = 4.58 x IO^^Mq 
and Rd = 4.0 kpc. Then the time unit is given by 24.9 Myr. 
We use a total of 6 x lO'' particles to simulate RCG004-A, 
RCG004-B and RCG004-C whereas for RCG004-0 a total 
of 10 X 10'' particles were used. The softening lengths for 
the disc, bulge and halo particles use d are 12, 5 and 20 pc 
respec tively following the suggestion of lMcMillan fc DehnenI 
(120071). These simulations are performed using the Gadget 
code (jSpringel et al.ll200il ) with a tolerance parameter 6toi = 
0.7, integration time step ~ 0.4 Myr. All the models are 
evolved for a time period of ~ 3.0 Gyr. 



3 SIZE AND GROWTH OF BARS AND 
BOXY/PEANUT BULGES 

Non-linear growth of a bar depends on the amount of an- 
gular momentum transferred from the inner region of the 
stellar disc to the outer disc, to a classical bulge (if exists) 
and to the surrounding dark matter halo in the host galaxy, 
all of which act like a sink. In this section, we focus on the 
effect of a rotating ICB on the size and growth of a bar and 
the final BP bulge. 

All 4 models described in section[2]rapidly form a bar in 
the disc. In the upper panel of Fig. O we show the evolution 
of the bar amplitude (defined as the maximum of the m = 2 
Fourier coefBcient A2 normalized to the axisymmetric com- 
ponent Ao) as a function of time for the four models. The 
net amplitude of the bar is lower in models with classical 
bulges having higher (Vm/a-)* . Such a decrement in the bar 
amplitude correlates well with the transfer of angular mo- 
menta to the bulge as shown in the bottom panel of Fig. [2] It 
shows the time evolution of the net gain of specific angular 



momentum normalized to the disc's initial angular momen- 
tum. The fast rotating ICB absorbs less angular momentum 
after ~ 1 Gyr as compared to the moderately rotating or 
non-rotating ICB. The most striking difference is between 
the fast rotating and non-rotating ICBs; after ^ 1 Gyr, the 
non-rotating ICB gains ~ 1.5 times more angular momen- 
tum than the fast rotating one (see Fig. [2]). Note that the 
linear growth rate of the bar is nearly unaffected by the 
initial rotation of the classical bulge. 

Fig. [3] shows the surface density maps of the inner 2Rd 
regions for all the 4 galaxy models at the end of 2.9 Gyr. 
In the face-on projection, they show strong bars and in the 
edge-on projection BP bulges. We measure the bar size using 
two methods: one using the drop in the ellipticity profile 
obtained using the IRAF ellipse fitting task on the FITS 
image (see Fig. [3| generated from the particle model. In 
the second, we use the radius at which the correspondin g 
position angle of the bar deviates by ~ 5deg (|Erwinll2005l ). 
We quote the bar sizes (denoted by Rtar) as the minimum 
of these two measurements obtained at the end of 2.9 Gyrs 
(se e Table[T]). For a list of methods on bar size measurement, 
see lAthanassoula fc Misiriotij (|2002l ). We have also used an 
independent method to derive the bar size directly from the 
particle model. In that, we use the FWHM of the radial 
profile of the bar amplitude as one measurement and the 
radius at which the phase angle of the bar deviates by 5 deg 
as the second. The minimum of these measurements is nearly 
the same as that obtained from the ellipse fitting method. 
From these measurements, it is found that the bar sizes are, 
in general, smaller if the ICB is rotating. On average, for the 
low mass ICBs considered here rotating with a (Vm/o')* ~ 
0.42 — 0.72, we found about 10 — 40% reduction in the bar 
size compared to the non-rotating case. 

The size of a BP bulge is calculated by finding zeros of 
the function Dg{x,z) defined as 

T,ios{x,z) -'Eios{0,z) 

^«("'^) = EU^) — ' « 

for a set of smoothed surface density (Ejos) profiles (see 
Fig. [4]) along slits parallel to the major axis (x) of a galaxy 
model such as shown in the lower panels of Fig. O 

The shape of the function Dg varies along the vertical 
direction. Close to the midplane (i.e., z = 0) of the galaxy, 
Dg is convex with the zeros occurring at a; = (i.e., along the 
minor axis) and values of Dg are always negative at a: 7^ 0. 
Slightly above the midplane (e.g., z = 0.11i?d in Fig. [Jjl 
Dg = in the inner regions {x ~ OARd) of some models 
and smoothly goes to negative beyond that. We define such 
a profile as characteristic of a boxy shaped bulge. At higher 
values of z (e.g., ai z = 0.23,0.3i?d, in Fig. |4| the shape of 
the function Dg takes the form of a characteristic peanut 
shape. Above z = 0.3Rd, we find that shape of the function 
Dg becomes irregular for some models and makes it harder 
to compare them with other models. We consider the profiles 
at z = 0.3Rd to compute the length of the BP bulge for 
all the m odels. The len g th of the BP bulge (denoted by 
BPL, see iLiitticke etHl (|2000l )) is obtained by calculating 
the distance between two zeros of the function Dg on either 
side of the galaxy centre and they are enlisted in Table [T] 
This method reliably returns the values of BPL as a function 
of vertical height (z) when there is a clear peanut signature 
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Figure 3. Bars and boxy/pcanut bulges at the end of 2.9 Gyrs. From left to right, the models are RCG004-0, RCG004-A, RCG004-B 
and RCG004-C. Initial rotation of the elassieal bulge inereases from left to right. Bottom panels show side-on projeetions of the models, 
above, along the short axis in the galaetie plane. 



and it is found that the largest BP bulge forms in model 
RCG004-0 with non-rotating ICB. 

The strength of the BP bulge in our models is obtained 
by calculating the maximum of the m = 2 Fourier coefficient 
(denoted by CJI""^) of the z-coordinates of the disc particles 
belon ging to the bar (|Martmez-VaIpuesta fc Athanassoulal 
|2008| V In order to understand the effect of ICB on the disc, 
we have used only disc particles to compute A2 and C2,z- 
The values of C™^ are shown in Table[T] It is found that the 
strength of BP bulges in our models with non-rotating ICB 
and fast rotating ICB are nearly the same. However, with 
moderately or slowly rotating ICBs, the strength is lowered 
by ~ 20 - 40%. 

Further, we perform a 2D bulge-disc decomposition on 
the F ITS images produ ced from our simulations using GAL- 
FIT (|Peng et al.ll2002l ) and derive the bulge half mass radii 
(^6,1/2), sersic indices (n;,) (see Table[T} and ellipticities for 
all the four final BP bulges shown in Fig. (3] Note that the 
BP bulge here is a composite bulge which is a superposition 
of the BP bulge formed out of the disk material and the ICB 
at 2.9 Gyrs. The initial value of Rb,i/2 for the non-rotating 
ICB is 0.26Rd- We do not find any definite trend either in 
the sersic indices or bulge sizes with the initial rotation in 
the classical bulges but the slowly rotating one in our model 
(RCG004-A) forms shorter bulge and have larger sersic in- 
dex. 



4 KINEMATICS OF BOXY/PEANUT BULGES 

Kinematically BP bulges are distinguished from the clas- 
sical bulges in that BP bulges possess cylindrical ro- 
tatio n as observed in a numbe r of external galax- 
ies ([Kormendv fc lUingworthl 1 19821 : iFalcon-Barroso et al.l 
l2006bl ) as well as in our Milk y Way (|Howard et all 120091 ') 
but see IWilliams et al.l |20lJ) for exceptional cases which 
displayed noticeable deviation from cylindrical rotation. All 
four models of BP bulges, studied here, rotate cylindrically 



Table 1. Rotating ICBs and properties of bars and boxy /peanut 
bulges at 2.9 Gyr in our simulations. Symbols are explained in 
seetion 2 and 3. 



Galaxy 




^bar 
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(xRd) 




(xRd) 




RCG004-0 


0.0 


1.35 


0.32 


1.12 


1.02 


0.085 


RCG004-A 


0.42 


0.95 


0.25 


1.33 


0.74 


0.056 


RCG004-B 


0.52 


1.20 


0.27 


1.13 


0.83 


0.068 


RCG004-C 


0.72 


1.10 


0.33 


1.03 


0.86 


0.084 
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-0.6 -0.3 



0.3 -0.6 -0.3 
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Figure 4. Shape of the funetion Dg along slits parallel to the ma- 
jor axis for all the four models. In each panel, profiles are drawn 
at z/Rd =0.009 (solid blaek line), 0.11 (green dashed-dot line), 
0.23 (blue dotted line), and 0.3 (red dashed line). The profiles at 
z/Fi.d = 0.3 are used here to measure the size of the boxy/peanut 
bulges. 
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but the degree of cylindrical rotation is seen to be vary- 
ing with time as the model galaxy evolves. In section 14.11 
we provide a simple prescription to quantify the degree of 
cylindrical rotation. The current section is focused on un- 
derstanding how the kinematics of the BP bulges depend on 
the initial rotation in the ICBs. 

In Fig. O we show mean LOS velocity profiles along 
slits placed parallel to the major axis of our galaxy model 
to probe how the rotation of a BP bulge varies along the 
vertical direction. At the end of 2.9 Gyrs, all the models ex- 
hibit clear (see section \AA\ for numbers) cylindrical rotation 
up to a height of ~ 0.23-Rd (~ Rb,i/2) which is about 900 
pc for the adopted scaling above the galactic midplane. Un- 
like models with rotating ICBs, the model with non-rotating 
ICB (RCG004-0) show cylindrical rotation up to a height of 
~ 1200 pc above the midplane. This is also clear from the 
2D velocity contours of the BP bulge region as depicted 
in Fig. [51 In models with rotating ICBs, the velocity con- 
tours start deviating above a height of ~ 0.257?d. Fig. [6] 
reveals that, although the initial bulge rotation for all the 
four models are quite different, at the end of 2.9 Gyrs they 
show similar kinematics. It remains to verify whether their 
evolutionary paths are also nearly identical. We address this 
issue in section [JT] 

In Fig. [71 we show the variation of mean velocity dis- 
persion (aioa) of stars in the final BP bulge along the ver- 
tical direction at two different radii from the centre of the 
respective galaxy model. In all the models, the minor axis 
c^ios profiles show steeper vertical gradient than the profiles 
a.t X = 0.17Rd (which is well within the bulge half-mass 
radii). A vertical gradient in the stellar velocity dispersion 
is observed in classic al bulges as well as in boxy bulges 
(j Williams et al.ll201ll ). It has also been observed that the 
velocity dispersion is remarkably constant along vertical di- 
rection in NGC4594 (a cl assical bulge) and in the well kn own 
boxy bulge of NGC4565 l|Kormendv fc Illingworthlll982l ). So 
a gradient in the velocity dispersion does not necessarily in- 
dicate the nature of the underlying bulge. If measured in the 
boxy bulge region, one would perhaps expect a shallow gra- 
dient in the velocity dispersion. From Fig. [71 it is clear that 
the rotation of the low mass ICBs does make little difference 
in the velocity dispersion of the final BP bulge. 

4.1 Criteria for cylindrical rotation 

By definition, cylindrical rotation means that the mean 
LOS velocity Vios{x,z) of stars at a particular distance 
(x) from the centre is independent of its vertical height 
(z) i.e., zero vertical gradient {dVios{x, z)/dz — 0). This 
is approximately true for observed ba rred galaxies as well 
as for bulges in N-body simula tions (|Combes et al.l 119901 : 
lAthanassoula fc Misiriotiil2002l '). Here, we construct a sim- 
ple criterion for measuring the degree of cylindrical rotation 
in BP bulges. 

We consider again a set of slits, nearly covering the ver- 
tical extent of the BP bulge, parallel to the major axis in 
edge-on projection. Along these slits, we have the LOS sur- 
face density T,ios{x,z) and velocity Vioa{x,z) of stars. For 
a given slit placed a.t z = Zi, we can compute the mass 
weighted line-of-sight velocity for all the stars within a pro- 
jected distance Xj from the minor axis of the edge-on im- 
age. This quantity may be normalized by replacing the in- 




-100 



Figure 5. The mean LOS velocity profiles along 4 different slits 
placed at z/R^ =0.009 (solid black line), 0.11 (green dashed-dot 
line), 0.23 (blue dotted line), and 0.35 (red dashed line) are shown 
in each panel. All the velocity profiles are shown for T = 2.9 Gyr. 
Cylindrical rotation is evident within the bulge half-mass radii of 
all models and extends beyond Ri,^i/2 the model with non- 
rotating ICB (RCG004-0). 




Figure 6. Spider diagrams for the mean LOS velocity of stars in 
the boxy bulge region which includes stars from the disk as well. 
From top to bottom, panels are drawn for models RCG004-0, 
RCG004-A, RCG004-B and RCG004-C. From left to right pan- 
els, they show evolution of the velocity structure between initial 
(T=0) and final states (T=2.9 Gyr). 
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Figure 7. The vertical profiles of LOS velocity dispersion at dif- 
ferent locations from the centre of the galaxy. The upper panel 
shows the dispersion profiles along the minor axis while the bot- 
tom panel shows for X = 0.17/?^- In both panels, solid black lines 
denote RCG004-0, green dash-dot lines denote RCG004-A, blue 
dotted lines RCG004-B and red dashed lines RCG004-C. All the 
dispersion profiles are derived at T = 2.9 Gyr. 



dividual velocities of all the stars within Xj with a veloc- 
ity at the bulge half-mass radius close to the disc midplane 
Vios{Rb,i/2, z ~ 0). After some experimenting, we adopt a 
weighting scheme and define a dimensionless quantity for 
the given slit as follows: 



5cL{z,Xj) 



ViosiRb,i/2,z ~ 0) /q^" T,ios{x,z)x^dx 



(2) 



We compute this quantity for the set of slits considered. 
The adopted weighting scheme puts more weight on the 
data points further from the minor axis. Note that, both the 
numerator and denominator have the dimension of angular 
momentum. If the stars close to the disc midplane rotate 
faster than the stars above, the denominator will always be 
greater than the numerator. If we had the knowledge of de- 
projection, one could interpret the above formulae as how 
the average rotational motion of stars associated with an in- 
finitesimally thin disc at a vertical distance z = Zi compares 
to the one close to the midplane. In general, the values of 
ScL are ~ 1 near the midplane and decrease as z increases. 

To proceed further, we plot the computed Scl values 
against the slit location {z) above the galactic midplane. 
The velocities of stars decrease along the vertical direction 
in our model galaxies and the same is expected for galaxies 
in equilibrium. So we expect Scl values can be well approx- 
imated by a linear function of z, where z is scaled to Rb.i/2- 
We then fit a straight line to the {Scl , z) curve and derive 



1 - 



0.5 



-0.5 



- 


' 1 




1 








1 

- 




e 


■ 




e 












e 




■ 










□ 










I 

▲ 


• 














★ 


■ 


• 






















□ 








★ 


■ 


ffi 


NGC 7332 

















NGC 4570 




□ 


















□ 






★ - 


• 


RCG004-0 








□ 








— RCG004-A 










□ 






RCG004-B 














■ 


— RCG004-C : 
1 


□ - 


-RCG004- 
1 


-C (T=0) 






□ 
1 



0.5 



1.5 



X/R 



b,l/2 



Figure 8. Degree of cylindrical rotation at different major axis 
radii of our N-body models of BP bulges plotted the end of 2.9 
Gyrs. Ovcrplotted here are the degree of cylindrical rotation for 
two galaxies NGC7332 and NGC4570. 
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Figure 9. Time evolution of the degree of cylindrical rotation 
calculated at the bulge half mass radii for four N-body models. 



the slope mcL- The degree of cylindrical rotation in the BP 
bulge is then defined as 



XXj) = 1 + mcz,. 



(3) 



The values of Scl bliq generally less than 1 unless Scl in- 
creases in the vertical direction which is unlikely for realistic 
models of galaxies. Note 52- z, is undefined for a non-rotating 
bulge because Eq. [2] would give ^ . Scl = 1 denotes pure 
cylindrical rotation. The degree of cylindrical rotation de- 
creases as the value of Scl decreases from unity. One could 
carry out the above exercise for different values of Xj (e.g.. 



0.5i?, 



Rb,l/2; '2Ri 



6,1/2 , -ft6,l/2 



'6,1/2 



etc.) and show how the de- 



gree of cylindrical rotation in a BP bulge would depend on 
the chosen value of Xj . 

In Fig. |8] we show how the degree of cylindrical rota- 
tion varies within the BP bulge at the end of T = 2.9 Gyrs. 
Within Rb,i/2, all the model BP bulges show nearly the 
same degree of cylindrical rotation. Beyond the bulge effec- 
tive radii, Scl decreases rapidly. For the models with mod- 
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erately rotating bulges (e.g., RCG004-A and RCG004-B), 
values are systematically lower (see Fig. [8|. We have 
also plotted the values of Sqj^ for the model RCG004-C at 
T = when the bulge showed no sign of cylindrical rotation 
(see Fig [6]). 

In order to test the above formulae, we applied it to two 
real galaxies, NGC7332 and NGC4570. NGC7332 is a well 
studied galaxy i n the literature and shows prominent cylin- 
drical rotation (|Falc6n-Barroso et al.|[20o3 ) . We have used 
SAU RON integral-field kinematics data (|Emsellem et al.l 
I2OO4I ') for the above two galaxies. The bulge effective radii 
{Rb,i/2) for NGC4570 a nd NGC7332 were taken to be 14 
and 6 arcsec respectively (|Emsellem et al.ll2004l : lFisher et al] 
I1994I ). Our analysis indicates that NGC7332 has a high de- 
gree of cylindrical rotation((52;^ ^ 0.8) within the bulge ef- 
fective radius. On the other hand, NGC4570 shows non- 
cylindrical rotation within its bulge effective radius (see 
Fig.©. 

Based on the above measurements of the degree of cylin- 
drical rotation in about 6 bulges, we introduce Scrit ~ 0.75 
as the boundary between cylindrical and non-cylindrically 
rotating bulges. Note that such a demarcation is not based 
on any rigorous calculation but should be treated as an op- 
erational definition. Scrit = 0.75 separates well, the rotating 
classical bulge (RCG004-C (T=0)) and non-cylindrically ro- 
tating spheroidal bulge of NGC 4570 from the rest of our 
sample. So a model bulge with Scl > ^cru {Scl < Scru) 
would be described as a cylindrical rotator (otherwise). It 
would be interesting to investigate 5qi^{X = Rb.1/2) for a 
larger sample of boxy bulges. 
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Figure 10. Top panels showing the surface densities for two 
models RCG004-0 (left) and RCG004-A (right) at T=1.05 Gyr. 
The middle panels show the contours of the 2D mean LOS ve- 
locity field for both the models. Velocity profiles at 4 different 
slits placed at z/Rd = 0.0(black solid line), 0.11 (green dash-dot 
line), 0.23 (blue dotted line) and 0.35 (red dashed line) parallel 
to the major axis of the galaxy are depicted in the bottom pan- 
els. Both models show deviations from pure cylindrical rotation 
within their respective bulge half-mass radii. 



4.2 Time dependent cylindrical rotation 

In Fig. [6l we showed that the kinematics of BP bulges at the 
end of 2.9 Gyrs turned out to be similar and independent of 
the initial rotation of the ICBs. Here, we compute the degree 
of cylindrical rotation within the bulge effective radius as a 
function of time for all the four models and plotted them in 
Fig. [21 This figure reveals interesting evolutionary behaviour 
for the BP bulges. First of all. Fig. [9] shows clearly that the 
degree of cylindrical rotation varies with time and that there 
are substantial differences at early times depending on the 
initial rotation of the ICBs. Nevertheless, irrespective of the 
initial condition, each model BP bulge evolves towards at- 
taining a higher degree of cylindrical rotation and within 
their respective bulge effective radii, the values of 5ql ^P- 
proach the same level. From Fig. [Ul we notice that it takes 
about 2.0 Gyr since the buckling episode for the cylindrical 
rotation to be present at the level of Scl 0.75 in our model 
bulges. In the section below, we show a snapshot during the 
evolution at T = 1.05 Gyr when cylindrical rotation in the 
bulge is still developing. 

4.3 Boxy/peanut bulges with non-cylindrical 
rotation 

Non-cylindrically rotating boxy bulges are seen in 
some galaxies e.g., NGC1381, NGC5746 and IC4767 
(|Williams et all 1201 ll ). It is not clear what might have 
caused such boxy bulges to rotate non-cylindrically. 

The analysis of the previous section demonstrate that 



the kinematics of BP bulges are time-dependent, especially 
the degree of cylindrical rotation. From Fig.[9l we see that at 
T = 1.05 Gyr beyond the buckling instability phase, two of 
the model bulges have less cylindrical rotation than the other 
two. In Fig. [To] and Fig. llll we show the surface density and 
kinematics of all four model galaxies at T = 1.05 Gyr. Ac- 
cording to our convention for Scrit, all 4 models of BP bulges 
at this time are non-cylindrically rotating. Note however, 
that the BP bulge in RCG004-C (see Fig. [11]) has cylindri- 
cal rotation within a smaller region about the minor axis. 
Since S^^iX — Rh,\/2) refers to the bulk of the bulge, such 
small details are averaged out. The present analysis indicates 
that boxy bulges with non-cylindrical rotation could still be 
in their early phase of evolution i.e., have formed relatively 
recently in their host galaxies. Another possibility is that in 
those galaxies the ICB was too massive to be spun-up all the 
way to cylindrical rotation; the response of massive ICBs to 
a growing bar will be investigated in a future publication. 

As shown in section [3] all the rotating and non-rotating 
ICBs gain angular momentum emitted by the bar. As a re- 
sult of the angular momentum gain, the ICBs spin up and 
their orbital structures evolve over time. A spherical har- 
monics analysis of the 3D density distribution of the ICBs 
showed that their inner regions formed a bar-like structure. 
For details on the case of non-rotating ICB (here , model 
RCG004-0), the readers are referred to ISaha et al] (|2012l ). 
The inner regions of these ICBs exhibit cylindrical rotation 
beyond about 0.5 Gyr which marks the bar buckling insta- 
bility phase in the disc. The net cylindrical rotation in the 
composite bulges (e.g., shown in Fig. [3] Fig. [10] or Fig. Ill|) 
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Figure 11. Same as in Fig. [TO] but for models RCG004-B (left) 
and RCG004-C (right). 



has its contribution from the evolved rotating ICB and the 
boxy bulge that form out of the bar buckling instability. 



5 CONCLUSIONS 

We have investigated the morphology and kinematics of a 
BP bulge that results from the instability of an axisymmetric 
disc assembled around a rotating low-mass ICB during the 
galaxy formation process. The composite BP bulge in our 
model contains a rotating ICB and shows highly evolving 
kinematics. Our main conclusions are highlighted below: 

1. The size and strength of the bars and resulting BP 
bulges in our simulations are reduced by ~ 10 — 40% in 
galaxy models containing ICBs rotating with (Vm/o")* = 
0.4 — 0.7 compared to one with (VWi/c)* = 0.0 . 

2. Rotating ICBs gain less angular momentum com- 
pared to the non-rotating ICB during the secular evolution. 

3. We show that cylindrical rotation in BP bulges in 
our simulations are time-dependent. During the early phases 
of BP bulge formation, significant deviations from cylindri- 
cal rotation can be observed. We quantify such deviations 
by providing a simple formulae. Our formulae when applied 
to NGC7332, gives a value for 5qj^ = 0.8 within the bulge 
half-mass radius. Using the same formulae, we confirm that 
NGC4570 is not cylindrically rotating. 

4. The degree of cylindrical rotation in the composite 
BP bulges does not depend on the angular momentum of 
the ICBs sufficiently long time after its formation. 

5. The early kinematics of BP bulges with rotating low- 
mass ICBs in our simulations suggest that the BP bulges 
with deviation from cylindrical rotation might have formed 
relatively recently. 
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